Modulation of Majorana induced current cross-correlations by quantum dots 
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We study charge transport through a topological superconductor with a pair of Majorana end 
states, coupled to leads via quantum dots with resonant levels. The non-locality of the Majorana 
bound states opens the possibility of crossed Andreev reflection with nonlocal shot noise, due to 
the injection of an electron into one end of the superconductor followed by the emission of a hole 
at the other end. In the space of energies of the two resonant quantum dot levels, we find a four 
peaked clover-like pattern for the strength of noise due to crossed Andreev reflection, distinct from 
the single ellipsoidal peak found in the absence of Majoranan bound states. 
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Majorana bound states (MBSs) are zero-energy 
fermionic states which are their own anti-particles. Since 
quasi-particles (QPs) in superconductors (SCs) are al- 
ways superpositions of electron and hole components, the 
Majorana criterion can be realized in a peculiar way: a 
zero energy QP in an SC has equal contributions from 
electrons and holes, and hence an exchange of particle 
and hole components leaves the QP invariant. There is 
currently much interest in the physics of MBSs [THE] , 
since one pair of MBSs nonlocally encodes a qubit, which 
is the building block for fault-tolerant topological quan- 
tum computing architectures. [f?l [TO] 

There is a variety of candidate systems for realizing 
Majorana fermions. Early proposals considered time- 
reversal symmetry broken p-wave SCs with the candidate 
Sr2RuO4.pl] Recently, the superconducting proximity 
effect has been suggested as a way to effectively induce 
p-wave pairing in topological insulators |12j and semi- 
conductors with strong Rashba spin-orbit coupling. |13L 
116] Recent experiments reported evidence of Majorana 
bound states (MBSs) in semiconductor-superconductor 
heterostructures.[TH5] A possible probe for the nonlocal 
nature of MBSs is crossed Andreev reflection (CAR), the 
conversion of an incoming electron into an outgoing hole 
in a different lead, [T?l - fS2] in contrast to local Andreev 
reflection (LAR), where electron and whole reside in the 
same lead. It has been shown theoretically that at suffi- 
ciently low voltages and small level width, LAR by a pair 
of MBSs is fully suppressed in favor of CAR.[23M] For 
voltages larger than the MBS energy splitting cm how- 
ever, local scattering processes dominate and suppress 
the crossed noise. 

Motivated by recent measurements of nonlocal trans- 
port in s-wave SCs coupled to quantum dots (QDs) 
[251 r25] . in this letter we study non-local transport 
through a coupled dot-MBS system as sketched in Fig.[T] 
In this setup, the QDs are used to suppress LAR due to 
Coulomb blockade. We consider a single-channel wire 
with one MBS located at each end. Each MBS is cou- 
pled to the adjacent lead via a quantum dot in the 



Coulomb blockade regime, modeled by a level of width 
F. Due to the finite wire length, the MBSs are tunnel 
coupled to each other and have an energy splitting ej\/ ~ 
Aexp(— L/£sm)> where £sm is the coherence length in 
the semiconductor. This tunnel coupling makes nonlocal 
CAR possible and gives rise to current cross-correlations 
oc (e 2 //i)e| / /r. In the regime T <C cm, transport occurs 
via discrete resonances, and CAR is suppressed if one dot 
level is resonant with the SC. This is due to the forma- 
tion of a zero-energy Majorana mode residing in the res- 
onant dot, |30| which cannot contribute to current cross- 
correlations since it has no energy splitting. The second 
regime r > ej; can be more easily realized experimen- 
tally. For large voltages eV cm, the crossed current 
correlator now has a four-leaf clover feature as a function 
of the dot levels ej, and e^. The extension of this struc- 
ture is determined by T, with a suppression of the noise 
when at least one dot is resonant with the SC. In strik- 
ing contrast, the crossed noise for two dots with a non- 
resonant proximity coupling shows a single-peak struc- 
ture with dominant weight along = — e_R.[Hl HH 
These findings are in excellent agreement with results for 
a microscopic model of a spinless p-wave superconductor. 

m 

Model system. — We consider the Hamiltonians 

H D = (^i d \ d i + 9i4^i + 9*Md^ (la) 

i=L,R 

HM = t M ilLlR+ ^ (Udl~fi+t*"fidi), (lb) 

i=L,R 

H s =&{d\d\ + d R d L ). (lc) 

Here, Hp describes two QDs coupled to leads, where di 
annihilates an electron with energy ej on dot i , tpi an- 
nihilates a lead electron, and gi is the lead-dot coupling 
strength. The lead electrons are characterized by their 
density of states pi , which is assumed to be energy inde- 
pendent, and have a chemical potential eV. We consider 
the regime where the QD single particle level spacing 
Se > eV > k B T, assume that the spin degeneracy is 
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FIG. 1. (color online) Schematic setup for a system with a pair 
of Majorana bound states (red dots) coupled to quantum dots 
which themselves are coupled to lead electrodes. The leads 
are biased with the positive chemical potential eV. Crossed 
Andreev reflection can be detected by correlating the currents 
Il and Ir that flow into the SC nanowire via MBSs. The 
nearby s-wave SC also induces a proximity pairing A between 
the dots. 



lifted by an external magnetic field, and that the QD 
ground state has an even number of electrons. Then, 
Kondo physics is absent and in the Coulomb blockade 
regime inclusion of only a single dot level in Hp is jus- 
tified. Hm describes two MBSs with an energy splitting 
cm coupled to the dots. The MBSs are described by her- 
mitian operators % — 7.J, which have anti-commutators 
{lit 7j} = 2<$ij, and are coupled to QD i with ampli- 
tude ti . The chemical potential of the superconducting 
wire hosting the MBS is zero. Hs describes an addi- 
tional proximity induced pairing between the dots with 
an amplitude A ~ 75 sin(fci?L)exp(— L/^sc)/{kFL) [21], 
where 7s is the normal-state QD level broadening due to 
the coupling between SC and QD, kp the Fermi momen- 
tum, L the wire length, and £sc the coherence length 
of the SC. We have in mind that this term may mainly 
be due to a coupling between the dots and the s-wave 
superconductor in a hybrid structrure. 

We diagonalize the Hamiltonian for MBSs and QDs 
without lead coupling by solving the corresponding BdG 



equation hM = e(l D + |Im)* with 



h = 



I ie M t L 
ie M t R 

t* T 
t 



' L 





1 R 



e L 
e R 






-A 



\ 

-t 
A 




'R 



-t L 



-e L 



(2) 



\ -t R A -e R J 



in the basis {j L , Jr, d\, dJ R , d L , d R }. Here, I D (I M ) de- 
note the identity matrix in the dot (Majorana) space. In 
the case A = 0, the quasiparticle energy spectrum has 
levels at 2cm, er> and e L , with avoided crossings where 
these levels intersect each other. If one of the QD energies 
is zero, the spectrum has a zero-energy MBS. 

To compute the zero-frequency noise through the 
above normal-SC-normal (NSN) system, we use a scat- 
tering matrix approach which also allows for Andreev 



reflection processes, 
noise correlators 



This yields the current and the 



ol k;y 



k,l;j,5 



(4) 



where Greek indices denote electron (e) and hole (h) 
channels, sign(e) = +1 and sign(/i) = —1, Latin indices 
denote the left (L) and right (R) lead, and 



,(ia) 



a/3* a?7 
S i,k S i,l ■ 



(5) 



The reservoir distribution functions Fermi func- 

tion with different chemical potentials for the electron 
and hole bands nk,-y = 1/(1 + exp(/3(e - sign(7)eV fe )). 



For the setup Fig. 1 
are the elements oft 

S(e) = 1 - 27nW f 



Vl = V R = V. The coefficients s 
the S'-matrix 



a,0 



-h 



w, 



(6) 

where W describes the coupling between the states of 
the system without leads and the scattering states in the 
leads, and [e(Io + ~~ h + inWW^]^ 1 is the retarded 
electron Green function for the closed system with self- 
energy inWW^ . The coupling matrix W in the lead basis 
WliVri ^l^r} is given by 



W iiai 



sign(a ld ) i gi I% /p 4l <5 iliid (5 cti 



(7) 



where ai d (a^ ) denotes the particle species of QD id (lead 
ii). The coupling strengths gi give rise to the level broad- 
ening Ti — 2irpi\gi\ 2 in the dots. In the following, we 
consider the case T L = T R = T ', t L = t R = t, and take 
the limit of zero temperature. 

Weak dot-lead coupling. — We begin our analysis in the 
regime A = and V < t < em . In Fig. [2j both differen- 
tial conductance and crossed current correlator Sl R are 
displayed as a function of bias voltage for several char- 
acteristic points in the e^-e^-plane. The differential con- 
ductance is peaked at the eigenenergies of Eq. (2). The 
peak width is determined by the broadening T. If one 
of the dot levels resides at the chemical potential of the 
SC, e.g. tL = 0, we always find one zero-energy state 
described by the Majorana operators 
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Here, 71 is localized on the resonant dot, while 72 is par- 
tially delocalized, and the weight of 72 on the resonant 
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FIG. 2. (color online) Current cross-correlator 5x.r in the 
weak dot-lead coupling regime with V — £m/4, t = 0.8eju, 
and A = 0. The lines for panel (b) are denned in (a), those 
for panel (c) in (d). The markers denote the results for the 
spinless SC model with e M = 0.01 meV and T = 0.002 meV. 



dot is determined by the energy cr of the non-resonant 
level. At the special point €l = £r — 0, we find 72 = 



/t 



In either case, the differential conductance has 



lR -ra R . 

a zero bias peak with height 4(e 2 / 'h)[t A / \t A + e 2 M e 2 R /4)]. 
Since the existence of a zero-energy MBS favors LAR over 
CAR, 23] we find that these resonances are dominated by 
LAR and yield only a small contribution to the crossed 
noise despite their large conductance. 

In contrast, we do not find a zero-bias conductance 
peak if both dots are non-resonant. In this regime, there 
is a striking difference between symmetric (e^ = €r) and 
anti-symmetric (ex, = —€r) positions of the dot levels. 
In both cases, we find contributions to the conductance 
and Slr due to the hybridization between the dots and 
the MBS. However, in the anti-symmetric case both the 
conductance and Slr are much larger than in the sym- 
metric case, and additional resonances at the QD energies 
contribute to crossed noise. This is due to the fact that 
Cooper pairs have zero energy, which leads to a suppres- 
sion of transmission through two resonant levels which 
have both the same energy in the symmetric case, but 
allows passage through QDs with opposite level energies 
in the anti-symmetric case. 

These findings agree very well with results for the 
microscopic model of a spinless p-wave SC defined in 
Eq. (10 1, sec Fig. [2j The only small deviation in Slr 
can be seen if both dots are resonant, where the effective 
model has a small negative Slr for large bias voltage, 
while it approaches zero for the microscopic model. This 
deviation has its origin in the presence of an additional 
transport channel due to a proximity coupling A in the 



microscopic model, which in principle could be described 
by the Hamiltonian H$ in Eq. (lc), but which is not in- 
cluded in the effective model H = Hm + Hr> considered 
here. 

Strong dot-lead coupling. — We consider the case t < 
Em <C T and begin with the situation A = 0. In 
Fig. [3^a), the correlator S L r for E M < eV = T/2 in 
the e^-e^-plane is shown. It is characterized by a four- 
leaf clover feature with a suppression of crossed noise 
along lines with either €l — or €r — 0, and peaks at 
= |cr| ~ T/2. While the peak height scales with 
e|f/r, the width of these peaks is larger than the Ma- 
jorana energy splitting due to the large value of T. As 
before, the suppression of the noise along t L = and 
cr = is mediated by the formation of zero-energy Ma- 
jorana modes by virtue of the dot-MBS coupling, which 
favor the LAR transport channel. 

The emergence of an approximate symmetry between 
symmetric and anti-symmetric positions of the dot levels 
(absent in the case T < t) can be understood as follows. 
For large T, the dots are strongly coupled to the leads 
and effectively become part of them. Hence, there are 
no separate resonances at the positions of the QD levels 
anymore, and only a single resonance due to the MBS in 
the wire survives. Since t <^L T, the broadening of this 
resonance is much smaller than T. As the QD levels can 
neither resolve this small broadening of the resonance, 
nor resolve the location of the resonance, the distinc- 
tion between symmetric and anti-symmetric QD levels 
becomes blurred, and the approximate symmetry arises. 
The Majorana zero energy state residing on one of the 
dots for £l — or cr — however does not change its 
character due to the presence of a large broadening T, 
and the noise stays low in this case, giving rise to the 
clover-like pattern in Fig. [3^a) . 

In Fig. (3^b), we complement these findings with re- 
sults for the microscopic model Eq. (10), for which a 



similar four-leaf clover structure emerges. However, sim- 
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FIG. 3. (color online) Current cross-correlator Slr for strong 
dot-lead coupling, (a) Effective model with eV = T/2, t = 
T/20, e M = T/10, and A = 0. (b) Spinless SC with e M = 0.01 
meV and T = 0.06 meV. For both (a) and (b), the pattern 
changes little for larger eV. 



4 




FIG. 4. (color online) Current cross-correlator Slr for eV = 
r/2. The dots are coupled (a) via the SC proximity effect 
with A = r/10 and (b) via SC proximity effect and with 
coupling to a pair of MBS with t = T/20, A = T/20, and 

£M = r/io. 

ilarly to the weak dot-lead coupling regime, there are 
small deviations with respect to the effective model near- 
er, = cr = 0, mediated by the SC proximity effect. 

To gain insight into the effect of an additional proxim- 
ity term Hs, we first discuss the situation without MBS, 
H = Ho + Hs- In Fig. |4|a), the crossed current corre- 
lator for the SC proximity case is shown. Here, Slr has 
a single peak of height oc A 2 /T near ej, = e# — 0, with 
width r along the direction cl = £r, and width eV along 
the direction cl — —£r- In contrast to the MBS case, 
there is no additional structure in this peak. 

In figure [4](b) , we consider the combined Hamiltonian 
H = Hm + Hd + Hs- We find a four- leaf clover feature 
similar to that in the Majorana only case, with the center 
of this feature now having a peak due to the proximity 
term in H$. From this, we conclude that the contribu- 
tions from the proximity effect and the MBS mediated 
CAR approximately add up. The relative peak heights 
in the crossed current correlator reflect the ratio of A 2 
and 

Microscopic model. — We complement our calculations 
by the analysis of a microscopic model for a spinless p- 
wave superconductor with Hamiltonian. [3T] 

JV-l N 

H K = - ^ (t K 4 +1 Cj+A K CjCj + i+h.c^) -ji K ^ c\cj, 

3 = 1 j 

(10) 

where the Cj annihilate a spinless fermion on site j 
with nearest neighbor hopping and nearest neighbor 
pairing amplitude A^. This model describes the low- 
energy physics of a nanowire in the topologically non- 
trivial phase. In the numerical analysis, we use the pa- 
rameters L = 1000 nm for the wire length, N = 200 
sites, tx = 20 meV, Ajf = 0.8 meV, and (Ik = 39.4 
meV, similar to the parameters used in |33j . These pa- 
rameter values yield the SC gap A sc = 0.3 meV and 
the Majorana energy splitting ejw = 0.01 meV. For the 
coupling of the operators ci and cat to the dots, we use 



tu.K = 0.025 meV. The results for this model agree very 
well with those for the effective model Eq. (1), see Figs. [2] 
and |3l 

Conclusion. — The nonlocality of a pair of Majorana 
bound states can be probed by crossed Andreev reflec- 
tion, whose observation is facilitated when suppressing 
local Andreev reflection with the help of two resonant 
QDs. In the case of a weak coupling between QDs and 
leads, we find a set of discrete transmission resonances. 
When at least one of the QD levels is tuned to the chem- 
ical potential of the superconductor, a zero energy Majo- 
rana state forms in the respective QD, which contributes 
only weakly to crossed Andreev reflection. This feature 
survives in the limit of strong dot-lead coupling, giving 
rise to a clover-like modulation of crossed shot noise as a 
function of QD energies, which is different from the single 
peak found without Majorana states. 

While finalizing this manuscript, we became aware of 
the preprint [34], where noise in a setup similar to ours 
is studied using a rate equation approach. 

We acknowledge helpful discussions with A. Das, 
M. Hciblum, and M. Horsdal, as well as financial sup- 
port from BMBF. 
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